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We study theoretically the spin-dependent transport in ferromagnet /normal- metal/ferromagnet 

f> |' double tunnel junctions by special attention to cotunneling in the Coulomb blockade region. The spin 

accumulation caused by cotunneling squeezes the Coulomb blockade region when the magnetizations 

in the ferromagnetic electrodes are antiparallel. Outside the squeezed Coulomb blockade region, 

' we propose a new anomalous region, where the sequential tunneling in one of the spin bands is 

suppressed by the Coulomb blockade and that in the other is not. In this region, the tunnel 

magnetoresistance oscillates as a function of bias voltage. The temperature dependences of the 

tunnel magnetoresistance and the magnitude of the spin accumulation are calculated. 
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D , The spin-dependent transport in magnetic materials, 
in particular, the tunnel magnetoresistance (TMR) in fer- 
romagnetic junctions, has attracted much interest 
Recent advances in nano-lithography and thin film pro- 
cessing make it possible to fabricate very small double 
I , tunnel junctions called single electron transistor (SET), 
where the electrostatic energy of excess electrons in the 
island has a significant effect on charge transport, i.e., the 
Coulomb blockade (CB) 0]. In the CB region, sequential 
tunneling, where tunneling events in each junction occur 
independently, is blocked at T = by the electrostatic 
energy, and superseded by a more complex tunneling pro- 
cess called "cotunneling" ||||. In ferromagnetic SETs 
' 0, it was recently pointed out that TMR oscillates in 
sequential tunneling regime || and TMR is enhanced in 
ON ■ the CB region by cotunneling j|. 

In this Letter, we study 

theoretically the spin-dependent single electron tunnel- 
ing in ferromagnet /normal-metal/ferromagnet (F/N/F) 
double tunnel junctions shown in Figs. |l|(a) and (c), 
and find a new anomalous region intrinsic to the spin- 
accumulated system. For the antiferromagnetic (A-) 
'■^J ■ alignment, where the magnetizations of the right and 
£h ! left electrodes are anti-parallel as shown in Fig. |l|(b), 
O ■ electrons with up(down) spin are easy (difficult) to tun- 
nel into the normal-metal, and difficult (easy) to tunnel 
out of it, because the densities of states of the left and 
J^j ' right electrodes are different between up and down spin 
bands. This imbalance among the tunnel currents causes 
the spin accumulation JTo| |l2[] , when the spin-relaxation 
time is sufficiently long in the normal-metal. In the CB 
region, where sequential tunneling is blocked at T = 0, 
cotunneling causes the spin accumulation. 

As shown in Fig. [l](b), the spin accumulation in- 
creases (decreases) the chemical potential for electrons 
with up(down) spin. This shift of the chemical potential 
decreases the energy for adding(extracting) an electron 
with down(up) spin to(from) the normal-metal island. 
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Therefore the critical voltage, above which electrons with 
down(up) spin can tunnel into(out of) the normal-metal 
island, is lowered and the CB region is squeezed as shown 
in Fig. |[ On the other hand, the critical voltage above 
which electrons with up(down) spin can tunnel into(out 
of) the normal-metal island is raised by spin accumu- 
lation. Consequently, a new anomalous region appears 
between the lowered and raised critical voltages. We call 
this region the half CB region, since the sequential tun- 
neling in one of the spin bands is suppressed and that 
in the other is not. Cotunncling is essential for appear- 
ance of the half CB region. If cotunncling is not taken 
into account, the spin accumulation does not occur in 
the CB region and the boundary of the CB region does 
not change, since there is no tunnel current. In the half 
CB region, we find an anomalous oscillation in TMR as 
a function of bias-voltage. The crossover between cotun- 
neling and sequential tunneling || is also found in the 
temperature dependence of the TMR and the spin accu- 
mulation. 

The TMR in the ferromagnet/ferromagnet (F/F) sin- 
gle junctions has been extensively studied , 
where the tunnel resistance depends on the relative orien- 
tation of the magnetizations between the electrodes, i.e., 
parallel or antiparallel. When the electrodes are made 
of the same ferromagnetic-metal, the tunnel resistance 
for the ferromagnetic (F-) alignment, where the magne- 
tizations of the left and right electrodes are parallel, is 
given by 1/Rf oc {T> 2 m +V^ n ), while for the A- alignment 
1/Ra oc 2V]\j'D m , where T>m and V m are the densities 
of states for the majority and minority spin bands at the 
Fermi level, respectively. The tunnel magnetoresistance 
is written as 

Ra-Rf _ 2P 2 

R F 1-P 2 ' () 

where P = (T>m — D m )/(T>M + T> m ) is the spin polariza- 
tion of the electrodes. The key point is that the difference 
between the products of the densities of states for each 
spin band causes the TMR in F/F junctions. 

We now turn to F/N/F double tunnel junctions shown 
in Fig. 0(a). It should be noted that the TMR does not 
exist in F/N single junctions, because there is no dif- 
ference between the densities of states for up and down 
spin bands in the normal- metal. For F/N/F junctions, 
however, the TMR exists, if the spin flip process in the 
normal-metal island can be neglected. The origin of the 
TMR in F/N/F junctions is quite different from that for 
F/F junctions. We assume that the right and left elec- 
trodes are made of the same ferromagnetic-metal such as 
Ni, Fe, and mixed-valent manganites like Lao.7Sro.3Mn03 
and the central electrode is made of a normal-metal such 
as Al with sufficiently long spin relaxation time. 

We first consider what happens if we neglect the 
Coulomb charging energy. For the A- alignment, in- 
troducing the chemical potential shift for up(down) 
spin electrons, /i|(/ij), the spin-resolved tunnel cur- 
rents through each junction at T = are given by 
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7i T oc T> N V M (\eV - Hi), hi oc V N V m (\eV - hi), / 2 f oc 
T>^V m {\eV + /if), and I 2 i oc V N T> M (\eV + /i;). Here, 
I? at is the density of states in the normal-metal for each 
spin band, and subscripts 1 and 2 represent the left and 
right junctions, respectively. 

The shifts, /if and /if, are determined by the stabil- 
ity condition, i.e., /if = I 2 -\ ancl III = ^24. • Then, we 
have Hi = —Hi = \P e ^ ' ■ The tunnel resistance for the 
A-alignment is written as Ra oc (1/T>m + 1/V m )/'WN. 
For the F- alignment, the spin-resolved currents are bal- 
anced without spin accumulation. Therefore, the chem- 
ical potential for each spin electrons does not shift and 
the tunnel resistance is given by the usual manner as 
R F oc l/V N {V M +V m ). The TMR for F/N/F junctions 
is obtained as 

Ra-Rf = (Pm - v m ? = p 2 (2] 

R F *D M V m 1 - P 2 ' U 

which is half of that for F/F junctions in Eq. (jl|). Equa- 
tion (||) has also been obtained by Brataas et al. [ p"5| . 
The polarisability P for Ni, Fe, and Lap^Sro^MnC^ 
are, respectively, 0.23, 0.40, (X6| and 0.83 In- 
troducing these P values into Eq. (^), we estimate the 
TMR for F/N/F double junctions made of Ni, Fe and 
La .7Sro. 3 Mn0 3 electrodes to be 5.6%, 19% and 220%, 
respectively. 

Let us next examine the single electron transistor with 
a capacitively coupled gate as shown in Fig. [l](c). The 
transistor is made of the same materials as the F/N/F 
double junctions discussed above, and its size is small 
enough to observe the CB. For simplicity, we assume that 
the insulating barriers for junctions 1 (left) and 2 (right) 
are the same; we subsequently set C\ = C 2 = C. 

The energy changes due to the forward tunneling of an 
electron with spin a through the junction 1 (n — > n + 1) 
and the junction 2 (n — * n — 1) are given by E^\ri) — 
E c (l + 2n)- ^eV g -\eV + Ha and ^(n) = E c (l - 
2n) + ^eV g — ^eV — /Vj respectively, where n is the 
number of excess electrons in the normal-metal island, V g 
the gate voltage, E c = e 2 /2C s , and C* s = 2C + C g . The 
energy changes due to the backward tunneling, E^ a ^ (n) 
and E^\n), are given by E[~\ri) — E^~\ri) + eV and 

E 2 V( n ) = E itH n ) + eV > respectively. 

The total current through the double junctions |l is 
expressed as 

00 

I = ^ Pn ^i^o-' (n) - I 2 <ria> (n)] , (3) 

n— — oc a,a' — ll 

where p n is the probability of charge state n and Ijaka 1 {n) 
represents the current, where electrons with spin a tun- 
nel into the central island through the jth junction 
and electrons with spin a' tunnel out of it through 
the fcth junction. The probability p n is determined 
by the condition for detailed balancing: p n [Tj + \n) + 
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r< +) (n)] - Pn+i[T{ \n + 1) + r£ -) (n + 1)], where 
rf } (») = E.(l/e 2 ^a)^ ) (n)/(exp[^(r l )/T] - 1) 
are the tunneling rates of n — > n ± 1 in the rth junction. 
Ijaka' (n) is obtained by the "golden rule" [H as, 



W (») = Jo „ p / de 1 d e2 de 3 de 4 /(ei) [1 - /(e a )] /(e 3 ) [1 - /(e 4 )] 

S (ei - e 2 + e 3 - e 4 + A£) , (4) 



e 2 - ei + (n) + *7 (+) e 3 - e 4 + W + *7 (-) 

where Jo = 3^ Rk — h/e 2 is the resistance 
quantum, Rj a is the tunnel resistance of the jth 
junction for electrons with spin er, AE is the en- 
ergy difference between initial and final states, and 
7 (±) = \T,.T,j9j^\n)coth[E^\n)/2T] repre- 
sents the decay rate of the charge state nil with 

9j° = 2§r~ Notc that in thc limit of lar S c tun " 

nel resistance (<? J(T — * 0) where the cotunneling current 
is negligible, Eq. (||) reduces to the sequential tunneling 
current in the orthodox theory Q . 

For the F-alignment, the tunnel resistances Rj„ are 
given by R\-\ = i?2| = Rm and Rii = R21 = R m , where 
Rm (oc \/T> mVn) and R m (oc \jV> m T>fi) are the tunnel 
resistances for electrons in the majority and minority spin 
bands of the ferromagnetic electrodes, respectively. The 
spin does not accumulate, since Ii^ji(n) = Iiiji(n) for 
/if = /ij = 0. Therefore, the CB region is the same as 
that for the usual metallic single electron transistors as 
shown in Fig. |l|(d). 

For the A- alignment, the tunnel resistances are given 
by R11 = R21 = Rm and R11 = i? 2 | = Rm- If = Mi = 
0, then Ii-\ji{n) ^ /^^(n), which gives rise to the spin 
accumulation and non-zero shifts fi a . The shifts of the 
chemical potential satisfy the condition that fif = — /Zj, 
because we assumed that the left and right ferromagnetic 
electrodes are the same and the density of states in the 
normal- metal is constant. We introduce the symbol 6 = 
Hi = —fi[, which is determined by the stability condition, 
E n Pn [IiVl (n) - h\j\ (n)] = 0. We carry out the 
numerical integration of Eq. (Q) and obtain 5. Here we 

take Rm = 2Rk {Rm = j^Rmj- 

At T = and for low V where sequential tunnel- 
ing is blocked by the CB, cotunneling causes the spin 
accumulation in the normal metal island, leading to 
the increase (decrease) of the chemical potential of up- 
spin (down-spin) electrons by 5 as shown in Fig. |l|(b). 
Therefore, the energy changes due to thc single elec- 
tron tunneling are given by E^\n) = Ef ] (n) ± 5 and 

E^in) = J3p°(n) =F 6, where J5, (±) (n) are the energy 
changes for 5 — 0. The CB region for the A- alignment is 

now determined by e[ (n), E2 (n) > 5 and squeezed 
as shown in Fig. ||. Note that for the F-alignment, the 
CB region is given by E[ ±] (n), E { 2 ±] (n) > 0. 

Outside the squeezed CB region, electrons with 
down(up) spin can tunnel into(out of) the normal-metal 
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island. However, those with up(down) spin cannot as far 

as Ei\n), E^in) > —8. As a consequence, we have 
a new spin-dependent CB region, where the sequential 
tunneling in one spin bands is suppressed and the other 
is not. This "half CB" region is given by the condi- 
tion that -S < E[ ±] (n), E^\n) < 8. In Fig. g, the 
boundaries of this "half CB" region for the A-alignment 
with Ni, Fe and Lao.7Sro.3Mn03 electrodes are plotted 
for n = 0. The solid line indicates the CB boundary for 
the F-alignment and the shade represents the half CB re- 
gion for the A-alignment with Lao.7Sr .3Mn03 electrodes 
in Fig. ^. One can see that the half CB region increases 
and the squeezed CB region, which is surrounded by the 
half CB region, decreases as P increases, because 8 in- 
creases with P. 

The tunnel magnetoresistance, (Ra ~ Rf)/Rf, and 
the shift of the chemical potential 8 at T = and V g = 
are plotted in Fig. |. At V ~ 0, the TMR for Ni, Fe 
and Lao.7Sr .3Mn03 electrodes are about 7.5%, 26% and 
310%, respectively. The values are 35 ~ 40% larger than 
those in Eq. (||). 

As the bias voltage V increases from toward the 
boundary of the squeezed CB region, the tunnel current 
for the A-alignment increases more rapidly than that for 
the F-alignment and the TMR decreases. At the bound- 
aries of the half CB region, Ra jumps because the cur- 
rent due to cotunncling decreases rapidly and that due 
to sequential tunneling starts to flow. The same jump 
in Rp appears at eV/2E c = 1, which is the boundary 
of the CB region for the F-alignment. Therefore, the 
TMR oscillates in the half CB region as shown in Fig. [|. 
We also find that cotunneling enhances the TMR around 
eV/2E c ~ 2.0, because it suppresses the V dependence 
of the total current for the A-alignment. 

We have also studied the system with large tunnel re- 
sistance, Rm = IORk, in order to see what happens 
when cotunneling is suppressed. In Fig. ||, TMR and 8 
for the electrodes of Lao.7Sro.3Mn03 and junctions with 
Rm = WRk are plotted by the dashed lines. The size of 
the half CB region is 79% of that for R M = 2R K . 

The temperature dependences of TMR and 8 have been 
calculated at V g — and eV/2E c = 0.1 in the squeezed 
CB region. As temperature T increases, the sequential 
tunneling, which is exponentially suppressed at low T, is 
recovered. The results are shown in Fig. ^. One can see 
clearly the crossover between cotunneling and sequential 
tunneling near T/E c = 0.1 @. 

In conclusion, we have studied the spin-dependent 
transport in F/N/F double junctions. In the antiferro- 
magnetically aligned SET, cotunneling brings about the 
spin accumulation at low V, where sequential tunneling 
is blocked at T = 0. The spin accumulation causes the 
squeezing of the CB region. Outside the squeezed CB 
region, we have found a new anomalous region, where 
the sequential tunneling in one of the spin bands is sup- 
pressed and that in the other is not. In this "half CB" 
region, the TMR oscillates as a function of V. The 
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crossover between cotunneling and sequential tunneling 
is found in the temperature dependences of TMR and S. 
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FIG. 1. (a) F/N/F double junction. Insulating barriers 
are shown by the shaded rectangles, (b) The densities of 
states for the antiferromagnetic-alignment are schematically 
shown. fj,f and /xj, denote the shifts of chemical potential for 
up- and down-spins, respectively, (c) F/N/F single electron 
transistor with a gate, (d) The Coulomb blockade region for 
the ferromagnetic-alignment with n = 0. 




FIG. 2. The boundary of the CB region for the ferromag- 
netic-alignment is indicated by the solid line. Boundaries of 
the half CB region for the antiferromagnetic-alignment with 
Ni, Fe, and Lao.7Sro.3Mn03 electrodes are plotted by dotted, 
dashed, and dot-dashed lines, respectively. The half CB re- 
gion for Lao.7Sro.3Mn03 electrodes is indicated by the shade. 
The tunnel resistance for the majority spin band is taken to 
be Rm = 2Rk for all systems. 
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FIG. 3. Left: The tunnel magnetoresistance for Ni and 
Fe electrodes are plotted against the bias voltage V in the 
upper panel. The chemical potential shift 8 is plotted in the 
lower panel. The shaded area represents the half CB region 
for Fe electrodes. The tunnel resistance for the majority spin 
electrons is Rm = 2Rk, the gate voltage is V g = and T — 0. 
Right: The same plot for Lao.7Sro.3Mn03 electrodes. The 
solid(dashed) lines indicate the results for Rm = 2Rk(WRk)- 
The half Coulomb blockade region for Rm = 2Rk is indicated 
by the shade. The gate voltage is V g — and T = 0. 




FIG. 4. Left: The temperature dependence of the tunnel 
magnetoresistance at V g = and eV/2E c — 0.1 for Ni and 
Fe electrodes is shown in the upper panel. The temperature 
dependence of the chemical potential shift 5 at V g — and 
eV/2E c — 0.1 for Ni and Fe electrodes is plotted in the lower 
panel. Right: The same plot for Lao.7Sro.3Mn03 electrodes. 
In both figures the tunnel resistance for the majority spin 
electrons is Rm = 2Rk- 
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